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On 

We present a method of post-Newtonian expansion to solve the homogeneous 
(2J[), Regge- Wheeler equation which describes gravitational waves on the Schwarzschild 



X 



spacetime. The advantage of our method is that it allows a systematic iterative 
analysis of the solution. Then we obtain the Regge- Wheeler function which is 
purely ingoing at the horizon in closed analytic form, with accuracy required to 
determine the gravitational wave luminosity to (post) 4 -Newtonian order (i.e., order 
v 8 beyond Newtonian) from a particle orbiting around a Schwarzschild black hole. 
Our result, valid in the small-mass limit of one body, gives an important guideline 
for the study of coalescing compact binaries. In particular, it provides basic for- 
mulas to analytically calculate detailed waveforms and luminosity, including the 
tail terms to (post) 3 -Newtonian order, which should be reproduced in any other 
post-Newtonian calculations. 



1. Introduction 



Among the possible sources of gravitational radiation, coalescing compact bi- 
naries have become one of the most promising candidates which may be detected by 
the near-future laser interferometric gravitational wave detectors such as LIGO [1] 
or VIRGO. Among various others, one very crucial reason for this is because of 
the simplicity of physics involved during the orbital evolution of such a system. 
Except for the very last moment of the coalescence, the dynamics is dominated 
by gravity alone and is almost non-relativistic. The energy carried away by grav- 
itational waves per orbital period is small enough and the orbit spirals in very 
slowly. Thus the gravitational radiation from the inspiraling binary has a distin- 

[3] 

guishable waveform which may be easily detected. Then by matched filtering, we 
may be able to determine the orbital parameters, spins of the constituent stars, 
the distance to the source, etc.. Astrophysical importance of such measurements 
is needless to say. Consequently much effort has been recently made to the study 
of the inspiral stage of coalescing compact binaries! 41 

However, based on calculations of the gravitational radiation from a particle 
in circular orbit around a non-rotating black hole, which is valid when one body 
of the binary has mass much smaller than the other and each has no spin, Cutler 
et al. showed that evaluation of the gravitational wave luminosity to a post- 
Newtonian order much higher than what can be presently achieved will be neces- 
sary to construct templates effective enough to determine the source parameters. 
Then the same problem was studied with much greater accuracy by Tagoshi and 
Nakamura' 61 in a semi-analytic way. In this very interesting work they calculated 
the coefficients of post-Newtonian expansion of the gravitational wave luminosity 
to (post) 4 - Newtonian order (i.e., 0(v 8 ) beyond Newtonian) and claimed that the 
accuracy to at least (post) 3 -Newtonian order is required for the construction of 
effective theoretical templates. In addition, they found logarithmic terms in the 
luminosity formula at and beyond (post) 3 -Newtonian order. 

It is then a matter of great importance to confirm those coefficients they calcu- 
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lated and to understand the origin of logarithmic terms. In this paper, we present 
a method to systematically calculate the gravitational radiation in the test particle 
limit and show that an analytic derivation of those coefficients is indeed possible. 
Our method is an extension of the one originally developed by Poisson, in which 
the expansion parameter is Muj where M is the mass of a Schwarzschild black hole 
(or a non-rotating compact star) and u is the frequency of gravitational waves, 
but in a way which is more suited for a systematic post-Newtonian expansion. In 
particular, we find it can more adequately deal with the logarithmic ambiguity 
one encounters when one matches light cones in the near zone and far zone in the 
post-Newtonian expansion. 

The paper is organized as follows. In §2, we formulate a new post-Newtonian 
approach to the gravitational radiation emitted by a particle orbiting around a 
non-rotating black hole, which is governed by the Regge- Wheeler equation. We 
then give an iterative method to solve for the Regge- Wheeler function, X™ £ (r), de- 
scribing a purely ingoing wave at the black-hole horizon, which plays an essential 
role in evaluating the gravitational radiation from a particle. In §3, as a prelimi- 
nary, we analyze general features of the iterative solution of X l Jl. We clarify the 
required accuracy of X™fc in order to evaluate the gravitational wave luminosity 
and waveforms to 0(v 8 ) beyond Newtonian. Namely, we find the calculation to 
O ((Mw) 2 ) is necessary only for £ < 3 to derive the luminosity formula to 0(v 8 ) : 
but additionally for £ = 4 to derive the waveform. In §4, we give the solution 
to 0(Muo) for arbitrary £ and to O ((Muj) 2 ) for £ < 4 in closed analytic form. 
Finally §5 is devoted to conclusion. Some useful mathematical formulas used in 
the text are given in Appendix. Throughout the paper we use geometrized units, 
c = G = l. 
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2. Formulation 



The situation we consider is the limit in which one of compact binary objects 
has a small mass relative to the other and each has no spin. In this case, the 
larger can be approximated by a non-rotating black hole (with mass M) and the 
smaller can be treated as a test particle (with mass fi <C M) oribiting around it. 
Of course, except for the case the larger body is actually a non-rotating black hole, 
approximating it by a black hole fails when the orbit approaches the event horizon. 
However, as long as one is interested in the regime in which one can neglect the 
radiation going into the black hole, this should be a valid approximation. In fact, 
we will see below that as long as we focus on the radiation outgoing to infinity, it 
is so until we reach 0(v 18 ) beyond Newtonian. 

It is known that the gravitational radiation from such a system is described by 
the Teukolsky equation,' 91 but it is also known that it can be transformed to the 
Regge- Wheeler equation,' 10 ' 111 which has a much simpler form. Here we consider 
the Regge- Wheeler equation. The inhomogeneous Regge- Wheeler equation takes 
the form, 

d2 2 (, 2M\ f£(£ + l) 6AA1 v 

d^2 +u ~\r) ~)\ m = m ' ( ' 

where S w £ m is the source term determined from the orbit of the particle and r* = 
r + 2Mln(r/2M — 1) is the tortoise coordinate. The explicit form of S u £ m for 
bound orbits can be found in Ref.12). 

The solution X u £ m which satisfies the outgoing-wave condition at infinity and 
ingoing-wave condition at horizon is then given by 

oo 
— oo 

(2.2) 

G(r, r') = 6{r - r')X out {r)X in {r') + 6{r' - r)X in (r)X out (r') , 

j^in d j^out -\rout d j^-in 

dr* dr* 
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where X m and X out are the solutions of the homogeneous Regge- Wheeler equation 
with the boundary conditions, 



X 



in 



X 



out 



Ce- iujr * for r* -> -oo, 

A in e -icor* + ^outgiwr* for r * _^ 

. . . . , (2-3) 

B m e -iur + B out e iur for r * _^ _ 



for r* — > oo, 



respectively. Note that from the constancy of the Wronskian for any two solutions, 
we have 

|^m|2 |^owt|2 |Cp |^go«t|2 |^gm|2 y 

W = 2iuA in = 2iuCB out . 



(2.4) 



Usually one normalizes X m so that C — 1, but we leave it arbitrary here for later 
convenience. 

Now as long as we are interested only in the radiation emitted to infinity, it is 
clear from Eq.(2.2) that we only need to know the form of X m . That is, we have 

Xuj£m{r) -> M u f m e , Auim = 2iuA in ' 

as r — > oo. Then the wave form at infinity is expressed as 

h + -^h x ~^J2J ^™ e_W((_r,) -^tm(0, <p) , (2.6) 

ira 

where _2^£m is the spherical harmonic function of spin weight s = —2. From the 
above waveform, it is straightforward to calculate the gravitational wave luminosity. 
For a circular orbit, the gravitational radiation spectrum is discrete: 

Ajim = $(w ~ mVL)Ai m , (2.7) 
where Q is the angular velocity. In this case, the gravitational wave luminosity is 
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riven by 



dE 
~dt 



16tt 



2, ,2 



im 



A ir . 



47T 



(2.8) 



In order to formulate a post-Newtonian expansion, we rewrite the homogeneous 
Regge- Wheeler equation as 

Xf = 0, (2.9) 

where z = ur, z* = z + e\n(z — e) and e = 2Mu, and hereafter we attach the 
eigenvalue index £ to X m for clarity. In the case of a circular orbit at radius 
r = ro, we have Qtq = (M/tq) 1 / 2 = v and uj = mfi, where v is the orbital velocity. 
Hence zq = ujtq ~ v and Muj = (M/ro)rou ~ f 3 . Further the source term Su£ m 
has support only for z < -2o- 112 '^ Hence the post-Newtonian expansion corresponds 
to solving (2.9) recursively with respect to e and then expanding the resulting Xf 1 
in powers of z as well as to extract out its incident amplitude Af 1 at infinity to 
required accuracy. 

Now, since z* contains e in itself, a naive expansion of (2.9) would result in 
a complicated equation which involves many e-dependent terms. Furthermore, 
such an expansion is not quite desirable at both boundaries since X 1 / 1 ~ e~ %z at 
horizon (z* — > — oo) and X] n ~ e ±lz at infinity (z* — -> oo). In order to take into 
account these asymptotic behaviors at both boundaries, one would have to take z* 
as the independent variable. However, this is not desirable either since the post- 
Newtonian expansion assumes an infinitesimally small e, but z* is not analytic at 
e — > (z* — > — oo as z — > e, while taking the limit e — > first yields 2;* — > 2; which 
never approaches —00). One possibility to circumvent this non-analytic behavior is 

f One may also transform X^ n to the corresponding ingoing-wave Teukolsky function first 
and then integrate the inhomogeneous Teukolsky equation. In this case, the source term 
has support only at z = z . [S1 



d z 

dz* 2 



+ 1 



+ 1) 



3e 

i 3 
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to separate out e~ tz dependence from the dependent variable Xf 1 in the beginning 
and to adopt z as the independent variable. Incidentally, this also takes care of 
the incoming part of the asymptotic behavior at infinity. In a sense, the incoming 
light-cone is adequately respected in this way. It turns out this idea works pretty 
well. Setting 

X in = e -ieln( Z -e) z ^ z ^ (2J0) 

we find (2.9) reduces to the form, 



£_ 2d_ f £(£+1) 
dz 2 z dz 



ee 



dz 



1 d , 

z 3 dz 



izjl 



(2.11) 



Surprisingly, e appears only as an overall factor on the right-hand side but nowhere 
else. Hence it readily allows us to apply the expansion in terms of e: 



(»), 



n=0 



for which we have 



(2.12) 



d 2 2d, 

dz 2 z dz 



= e 



dz 



6%) 



(2.13) 



Once & is obtained to desired accuracy, Xf 1 is readily obtained by inserting 
the resulting ^ into Eq.(2.10). The calculation of AJ 1 is then straightforward by 
gathering together all the coefficients of leading terms proportional to e~ %z / z from 
the asymptotic form of ^ at z — oo. In the rest of this section, we describe a 
method to carry out this procedure. 

Solving Eq.(2.13) recursively is in principle straightforward. Since the general 
homogeneous solution to the left-hand side of it is given by a linear combination 
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of the spherical Bessel functions je and ri£, one can immediately write the integral 
expression for Noting that jgn' e — = 1/z 2 , we have 



4*^ J 



(2.14) 



where the prime denotes d/dz and we have tentatively expressed the integrals in 
the indefinite form for convenience. If these indefinite integrals can be explicitly 
performed, the boundary condition for Xf 1 as given in Eq.(2.3) is easily imple- 
mented. For this purpose, we first rewrite Eq.(2.14) by performing integration by 
parts as 

= ~nt J dz (zh-^'^^e^t^' + je J zdz [j e ^n £ ] 



-n e jdzz^'4 n ^' + kj e i 



7 f i*fi( n ~V' - i'^ n ~ V ' 



+ k I dz [j e -> m] 



= n e Jdz{ + je) it 1] + " (j^' ' fet*) } 

- je J dz [j e -> ne] + z(jtri e - n e j' e )^ n ~ 1] '. 

(2.15) 

Then using recurrence formulas for the spherical Bessel functions, the above is 
further rewritten as 



An) _ lJn-1) • An-1) _ j(n-l) 



+ in e J dzz(ji$ n 1)1 - j'gif l) ^j-ije J dzz(n&f ^'-n'^f l) ) , 



(2.16) 
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where is an indefinite integral defined as 
Here Q stands for either ji or ri£. 

(n) 

Now let us consider the boundary conditions for Q . The boundary condition 
for Xf as z* -> -oo (z -> e) is that Xf = e- £ln (*-e)2& oc e"* 2 *. Hence s& is 

(n) 

regular ai z — e. Since e can be taken arbitrarily small, this implies that zQ 
must be no more singular than 0(z~ n ) at z — 0. In particular, this fixes the lowest 
order solution to be £^ = A^j^ where is an arbitrary normalization constant. 
For convenience, we set Ng — 1. Taking account of the behavior of the lowest 
order solution and the structure of the basic equation (2.13), we then infer that 
z£^ must be no more singular than ^ +1 ~ n at z — 0. Thus for n < 3, which is 
sufficient for the purpose of this paper, the boundary condition is that all the zQ 
are regular at z — 0. In passing, we note a symmetry of ^. From the form of 
Eq.(2.13) and the fact that ffl = ji and ji{—z) = (—l) e ji(z), we deduce that 

&\-z) = (-lY+ n $%) , (2.18) 



where the bar denotes complex conjugation. It is worthwhile to keep this symmetry 
in mind. 

Since j £ = 0{z e ) as z -> 0, we have X} n -> 0(e £+1 ) e -^, or Q = 0(e £+1 ) in 
the notation of Eq.(2.3). On the other hand, from the asymptotic behavior of jt 
at z = oo, we find A™ and A° ut are of order unity. Then from the first of Eq.(2.4), 
we obtain 

I a in | i/i out I 

Thus |A^ n | = until we go to 0(t 2t+2 ) or more. Since £ > 2, this implies 

that even though we are dealing with the Regge- Wheeler equation which describes 



9 



gravitational radiation from a non-rotating black hole, it also correctly describes a 
situation with a general non-rotating compact star, provided that the inner struc- 
ture of the star may be ignored and we focus on the outgoing radiation up to 
the amplitude of less than 0(e 6 ) or 0(v m ) beyond Newtonian. We note that if 
one considers the effect of radiation reaction to the orbital evolution, difference 
appears between the cases of a black hole and a star much earlier, at 0(v 8 ) beyond 

[131 

Newtonian, due to the radiation absorbed into the black-hole horizon. What we 
claim here is that the effect of the presence of the horizon on the outgoing radiation 
appears at 0(v ls ). 

The above fact also implies that we may choose the phase of X| n so that Af 1 and 
A™ 1 are complex conjugate to each other. Then since the Regge- Wheeler equation 
is real, it follows that Xf 1 can be made real up to 0(e 5 ). Provided we choose the 
phase of X™ in this way, Im for a given n is completely determined in terms 

of Re (t^p j for r < n — 1. To see this, let us decompose the real and imaginary 
parts of 

= ^ + 4 n) ■ (2-20) 

Inserting this expression into Eq.(2.10) and expanding the result with respect to e 
by assuming z ^> e, we find 

X in = e -ieM Z -e) z ^ + <f (l) + ^(D } + £ 2 (/ (2) + ^(2) } + ^ 

=z (j t + e/, (1) + e 2 (/f + gf ] ]nz- ^(M 2 ) + • ■ •) (2.21) 

+ iz (eigP - j, In z) + e 2 (gf ] + \ - In z) + ■ ■ . 

Hence we must have 

gP=jt\nz, gf ] ^--Ji + ffhnz, (2.22) 

Since the reality of X^ n is not manifest in our expansion scheme, this fact can 
be used as a useful method of checking calculations. To conclude this section, we 
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note the relation between the functions fj> n ^ and the conventional post-Newtonian 
expansion of Xf 1 : 

oo 

vin \ n v{n) . 

n=0 



3. Preliminary analysis 



(2.23) 



Prior to actually carrying out the integration of Eq.(2.16), in this section we 

(n) 

investigate general features of the asymptotic behavior of Q at z = and z = oo 
for n — 1 and 2. In addition, we comment on the behavior of Q ' at z — for 
£ < 4. Consequently, we clarify the condition to achieve 0(v 8 ) beyond Newtonian. 

First we consider the asymptotic behavior at z — 0. We return to Eq.(2.13) 

(n) 

and represent the right-hand side of it by S) . As we mentioned at the end of the 

(n) 

previous section, the imaginary parts of Q are subject to their real parts of lower 
n. Hence we may focus on the real parts. 

For n — 1, we have 

for £ = 2, (3,1) 
) for £ > 3. 

Hence from Eq.(2.13) or (2.14) and noting the regularity at z — 0, we obtain 

^' ^ ' 1 O^" 1 ) + «('>i( 2 ) for£>3, 

where c/jp is an arbitrary constant. Since its presence simply corresponds to a 
different normalization of X l e n , we set ocp = for simplicity. 
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For n = 2, we then have 



Re (Sf >) = - ( /<«» + !/«< - 1±4<») - I ( 2,/P' + ^> 



1 1 jO(/- 2 ) forC = 2,3, 

^ V z 2 \o(z £ ^) for£>4, 



(3.3) 



where gf } = jg\nz as given in Eq.(2.22). As may be expected, the terms in 
Eq.(3.3) arising from this imaginary part just cancel the j^(lnz) 2 /2 term of 
in Eq.(2.23). An easy way to check this is to insert it into the left-hand side of 
Eq.(2.13). Hence keeping in mind the regularity condition again, we find 



, 0(z £ ) + 0(z £ ) In z- iMlnz) 2 for £ = 2, 3, 
1 0(z £ - 2 ) + 0(z £ )lnz- \je(\nz) 2 for £ > 4. 



Inserting Eqs.(3.2) and (3.4) into Eq.(2.23), we find 

X l 2 n = z 3 [0(1) + eO(z) + e 2 {0(1) + 0(1) \nz} + ■■■], 

Xf = z 3 [0(z) + eO(l) + e 2 {0(z) + 0(z) \nz} + ■■■}, 

X} n = z 3 [0(z £ - 2 ) + eO(z £ - 3 ) + e 2 jo(/- 4 ) + 0(z £ ~ 2 ) In 2} + • • •] (£> 4). 

(3.5) 

(n) 

One sees that X^ behaves more regularly than what may be expected from the 
boundary condition we discussed in §2. As noted previously, in the case of a 
circular orbit, we have z = 0(v) and e = 0(v 3 ). Then from Eq.(3.5), we find that 
0(e 3 ) terms for £ = 2 and 4 may contribute to the waveform at 0(v 8 ) beyond 
Newtonian, if the leading terms of them behave as z" 1 relative to the Newtonian 
term. By explicitly calculating the series expansion formulas for £ = 2 and 4, we 
find it is the case for £ = 2 but not for £ = 4; the leading term of order z~ x relative 
to Newtonian cancels out for £ = 4. For completeness, we give the leading term of 
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X2 here: 




319 



z 2 + 0(z 4 ) + 0{z 4 ) In z. 



(3.6) 



6300 



Thus in order to evaluate the gravitational waveform and luminosity to 0(v ) 

(n) 

beyond Newtonian, it is necessary to know the series expansions of/; ; for£ < 10- 

(3) 

2n with n — 0, 1, 2 and the leading term of ■ Note, however, that calculations 
of these series expansion formulas are straightforward and easy to obtain. Note 
also that, as can be seen from Eq.(2.8), if we are interested in the luminosity alone, 
we only need to know for £ < 6 — n up to n = 2 and the leading term of ■ 

Recalling again the general formula for the luminosity, Eq.(2.8), and the sym- 
metry of ^ n) , Eq.(2.18), we see that the In z terms in the behaviors of X™ and 
X| n should give rise to lnt> terms in the luminosity at 0(v G ) and 0(v 8 ) but not 
at 0(v 7 ) (note that the \nz appearing in the above expressions is, mathematically 
speaking, actually ln|z|). This confirms the correctness of the result obtained by 
Tagoshi and Nakamuraf'at least qualitatively. Although it may be too early to 
conclude anything, the present consideration indicates that the Inv terms in the 
luminosity originate from some local curvature effect near the source (local in the 
spatial sense; it may well be non-local in the temporal sense), but is not due to the 
effect of curvature during propagation. This last point will be discussed further 
below. 

Now we turn to the asymptotic behavior at z = oo. For this purpose, let the 

(n) 

asymptotic form of f) ' be 



Then noting Eq.(2.22) and the equality e t€ln ( z e ) = e lz * e lz , the asymptotic form 



in) p (n) 



'je + Qi n l as 2 ^ oo (n = 1, 2). 



(3.7) 
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of Xf 1 is expressed as 



J^in ieln(z— e) 



\e~* (zhfe*) \l + e{pW + i(QW + In z)} 

+ ^{(pf)-g«m,) + i(Qf + P«in,)} 

l + e^-i^-ln*)} 



+ 



+ -e xzr [zh { Pe- iz )e" 2ieln ^ 



(3.8) 



+ e 2 {(pf ) + gf ) ln,)-i(gf -P £ (1) ln ,)} + •• 



where Yip and /i^ are the spherical Hankel functions of the first and second kinds, 
respectively, which are given by 



= k + int - H) m - , hf = u - int - 



(3.9) 

4, 4, 

Using the above asymptotic behavior of h^p and the incident amplitude Af 1 
can be readily extracted out: 

= i^ +1 e- i£ln£ [l + 6 {p, (1) + i (q { P + Inz) } (3-10) 
+ e 2 {(pf -g« In,} +i(gf + P« In,)} 

where note that 



+ 



u;r* = o> ( r + 2Mln 



r-2M 
2M 



z* — e In e, 



(3.11) 



from our definition of z*, which gives rise to the phase — ielne of Af 1 . Since we 
have required Xf 1 to be real, we must have A° ut = AJ 1 . This can be explicitly 
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checked by expanding in powers of e the factor e 2leln ( z e ) associated with the 
terms proportional to e lz in Eq.(3.8). 

An important point to be noted in the above expression for AT is that it 

in) (n) 

contains In ^-dependent terms. Since AT should be constant, and Q\ should 
contain appropriate In z-dependent terms which exactly cancel the In z-dependent 
terms in the formula (3.10). To be explicit, we must have 













— \az , 


p(2) 




+ q { p In 


Qf 


= #> 


— ' m 



(3.12) 



where pf^ and qjfi (n = 1,2) are constant. These relations, as well as Eq.(2.22) 
mentioned previously, can be used to check the result of calculations. 

In terms of pf^ and q^ n \ AT is finally expressed as 



A £ — 2 



l + e( P [ 1) +iq { e 1) )+^(p^ + iq ^) + ..]. (3.13) 



Note that the above form of AT implies that the tail radiation, which is due to 
the curvature scattering of waves, will contain ln-y terms as phase shifts in the 
waveform, but will not give rise to such terms in the luminosity formula. This 
supports our previous argument on the origin of the liav terms in the luminosity. 

Since e = 0(t> 3 ), the accuracy of AT to 0(e 2 ) is sufficient for calculations up to 
0(v 8 ) beyond Newtonian. Specifically, taking account of Eq.(3.5) and the remarks 
that follow it, we conclude that in order to calculate waveforms to 0(v 8 ), the 
accuracy of AT we have to achieve for a given £ is 0(e n ) where n = [^] . On the 
other hand, if we concentrate on the luminosity, the required order is n = [ 12 g 2 ^ ] . 
Thus, in order to derive the luminosity formula accurate to 0(v 8 ), we need to 
evaluate AT for only i = 2, 3 to 0(e 2 ) and for £ = 4 to 0(e), but to obtain the 
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waveform to 0(v 8 ) we additionally need A 1 ™ for £ = 4 to 0(e 2 ) and for £ = 5 to 
0(e). From the observational point of view, however, the post-Newtonian terms in 
the luminosity formula are far more important than those in the waveform because 
of their accumulative effect on the orbital evolution of compact binaries? 1 

4. The ingoing-wave Regge- Wheeler function to 0(e 2 ) 

In this section we derive closed analytic formulas of the ingoing-wave Regge- 
Wheeler function for arbitrary £ to 0(e) and those for £ = 2, 3 and 4 to 0(e 2 ). 

4.1. Calculation to 0(e) 

From Eq.(2.16) with ffl = jg, we have 



1 



;3 Z 



Although a bit tedious, the integrals A 1 } and I* 1 ] can be straightforwardly eval- 
uated by using integral formulas for the spherical Bessel functions. Some useful 
formulas are given in Appendix. In addition, the fact that ji and ri£ appear in 
an anti-symmetric manner in Eq.(4.1) helps us to save our computational effort 
considerably. We then find 

(i) _ (£-!)(£+ 3) . _ / I 2 -4 21-1 \ . 
^ ~2(£+ l)(2£+ l) 3e+1 \2£(2£+l) + £(£-!)) Ji ~ 1 

(4.2) 



+ z 2 (n e j - j e n )j + V ( — + — 1 —) 

^— ' \m m + 1 / 

m=l N ' 



Z (ngjm ~ 3 mm) 3 7 



+ ri£ (Ci 2z - 7 - ln 2z) - ^ Si 2^ + ijg \nz + op ji + , 

where Si x = Jq dt sin t/t and Ci x = — dt cos t/t are the sine and cosine integral 
functions, respectively, and atp and (3^p are integration constants to be determined 
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from the boundary condition. Note that terms like z 2 (nij m — jg.n r 
line are polynomials in 1/z which are known as the Lommel polynomials: 

z 2 {n £ j m -jin m ) = Ri- m - 1)m +s(z), 



in the second 



r=0 



(l-m-l-r)\T(l+\-r) (2 
r\(£-m-l- 2r)\ T(m + § + r) 



£-m-l-2r 



(4.3) 



As we discussed in §2, the boundary condition is that z^jp is regular at z — 0. 
This fixes @jp = 0. As for af \ since it simply represents the arbitrariness of the 
normalization of X| n , we set ocp = for convenience (in Poisson's notation) 81 our 
choice corresponds to his ct£ = —n/2). Note that this choice is in accordance with 
the choice we made in the previous general discussion, as seen by comparing the 
asymptotic behavior of tffl at z = and Eq.(3.2). 

Asymptotic behavior of £^ at z = oo is easily seen. Noting the fact j£ + \ ~ 



k-\ ~ (-l) £+n n 2n -£, etc., we find 



-31 



p { p + i In z ) j£ + (<?f ) - In z j n £ ; 



(4.4) 



where 



7T 

2 ' 



9i 



(1) 



*/>(£)+ i/>(£ + l) + 



(I -!)(* + 3) 
£(£+1) 



-In 2. 



(4.5) 



Here ip(£) is the digamma function: 



f-1 



m=l 

and 7 = 0.57721 • • • is the Euler constant. Note that Eq.(4.4) is perfectly consistent 
with general discussion given in §3. A™ to this order is readily found to be 

Af = -i e+1 e- iehl€ \l + e(-- 
2 - V 2 



• (1) 



+ 



(4.6) 



with qj, given by Eq.(4.5). 
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Before closing this subsection, we mention an important issue associated with 
the interpretation of Eq.(4.6). Comparing our result for 1 = 1 with Eq.(4.24) of 

rgj 

Poisson's paper, we see that his lne term (e = 2Mcu) at 0(e) is automatically 
absorbed into the phase of A™ in our result. An inspection of Poisson's analysis 
reveals that his In e term indeed originates from the phase of Xf 1 when matching 
r* and r; r* = r + 2Mln[(r - 2M)/2M], in the asymptotic behavior of Xf 1 ; 
e ±icur _ e ±i(z -elne) Thus, as discussed in the previous section, this lne term will 
never contribute to the amplitude of A™ but only to its phase at any higher order of 
the expansion. Then since a different choice of a constant scale in the logarithmic 
term of r* as 2M — > 2M' induces only a change of phase by iem(M/M'), the lne 
term itself has no physical significance except for its In ^-dependence. 

The above argument implies that one should not take the value of qp itself as 
physically significant, since it may also be absorbed into the re-definition of r* in 
the phase of X 1 / 1 . The physical significance lies in the relative differences of 
between different I. In passing, we note that Wiseman found a discrepancy in the 
tail term of the waveform between his result of a post-Newtonian analysis' 141 and 
Poisson's result, which corresponds to a discrepancy in the value of q£\ Wiseman 
argues that this discrepancy may be due to different boundary conditions in the 
two approaches; with and without the black-hole horizon. Because of the reason 
we mentioned above, in addition to general consideration on the presence of the 
horizon we gave in §2, we suspect Wiseman's argument is incorrect. 

4.2. Calculation to 0(e 2 ) 

Using recurrence relations for the spherical Bessel functions, one can reduce 
all the terms in £^ of Eq.(4.2) to a linear sum of the spherical Bessel functions, 
except for the terms proportional to Si 2,2, Ci2^ and \nz. Hence if the integrals 
involving these latter terms can be evaluated when inserted into Eq.(2.16), we may 

(2) 

obtain an analytic formula for Q . Performing integration by parts and utilizing 
the anti-symmetric appearance of jl and ri£, we then find this condition translates 
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to the condition that the following integrals are analytically calculable: 





where for notational convenience we have defined the functions C and S as 

z 2z 

C(z) := J dx (-2xj (x) 2 )= J dx C °* X x ~ l = Ci2^- 7 -ln2z, 

* " <«) 



/r sin x 
(— 2xno(x)jo(x)) = / = Si2z. 





We note that C(z) = 0(z 2 ) and S(z) = 0(z) so that B^(z) are perfectly regular 
at z — 0. Now using the fact nj^z) = — Jq{z) + l/z 2 and the above definitions of 
C(z) and S(z), we can rewrite -B^(z) (( = n,j) as 



z z 

B n (z) = j dx (C'(x)C(x) + S'(x)S(x)) +2 J ^C(x) 





2z 



1 /" dr 

= (^W 2 + ^) 2 ) +2C(2)ln22 + 2 / — (1 -cosx)lnx, 

Z J X 



Bjy 



i{z) = J dx (C'(x)S(x) - S'{x)C{x)) (4-9) 


2z 

— sin x In x 



2z 

/dx 
— (Si x cos x — Cix sin x) . 
x 



Hence the asymptotic behavior of at z = 00 can be obtained if the integrals in 
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the above expressions can be evaluated for z — > oo. In Appendix we show this can 
be done: 

2z 

/dx 
— sin a; In a; — > —717 , 
x 



2 ! —{I - cosx) lux -> (ln2z) 2 + ^- - 7 2 , (4.10) 
J x 12 


22 

/G?X 
— (Si x cos x — Ci x sin x) — > 0. 
x 



Using these results, we find 



B n ^^ 2 - l -{ 1 + \n2z)\ 
£j-~( 7 + ln2z), 



(4.11) 



as 2; — >■ 00. 



(2) 

Thus, it is possible to derive an analytic formula of Q for arbitrary £, though 
the manipulation will be quite involved. However, as discussed in the previous 
section, as long as we are concerned with the accuracy up to 0(v 8 ) beyond New- 
tonian, it is unnecessary to derive a general formula of Instead, at 0(e 2 ), 

(2) 

we only need the series expansion formulas of Q at z — for £ < 6 and the 

(2) 

asymptotic formulas of Q at z = 00 for £ < 4. The former can be easily obtained 
by expanding in powers of z the integrands in Eq.(2.16). Hence it is sufficient to 

(2) 

derive closed analytic formulas of Q up to £ — 4 for calculation of the luminosity 
and waveforms with accuracy to 0(t> 8 ). In what follows, we present the results of 
our calculation for £ = 2,3 and 4 in order. 
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(a) t = 2 

First consider the case of £ = 2. From Eq.(2.16), we have 

/ 2 (2) =Re(ef) 

=\i l) + 32 j dz(2zn 2 - n 3 )/ 2 (1) - n 2 J dz(2zj 2 - j 3 )f^ 
+ 32 J dzn 2 j 2 -n 2 J dzj\ - ^j 2 (hi z) 2 , 
9 { ? =Im(«f) 

=-,72 In z + j 2 / dz(2zn 2 - n 3 )j 2 In z - n 2 / dz(2zj 2 - j 3 )j 2 In z 



(4.12) 



- h J dz z (n 2 f^ f - /<V) +n 2 Jdzz (j 2 f^> - /< V) . 
Calculation of the above integrals is straightforward but requires care and pa- 

(2) 

tience. As discussed in §2, the imaginary part should be equal to the one 
given in Eq.(2.22). With the boundary condition that it be regular at z — 0, we 

(2) 

have explicitly evaluated the above integral expression for and confirmed our 
expectation. 

As for the real part f^, we obtain 

,(2) 389 113 1 r-t / \ 
f2 =-^«-^-n + T n + Bn{z)32-B 3 {z)n 2 

N / 5 5 3 107 \ 
+ C{Z) \Vz n2 ~ 3 H1 ~ ~z n0 ~ 2T0 32 ) 

(4.13) 

„, , / 5 . 5 . 3 . 107 \ 
~ S{Z) {6z j2 -3 Jl --z J0 + 2T0 n2 ) 



107 ■ i 1 • n ^2 , (2) ■ 
~2W J2 2 32 ( ' 2 J2 ' 



(2) 

where C, 5, _B n and Bj are the functions defined in Eqs.(4.7) and (4.8), and a 2 is 
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an arbitrary constant which we shall set to zero for convenience. For confirmation, 

(2) 

we have examined the above expression for with the help of the algebraic 
manipulation program, Mathematica. We have confirmed that it satisfies the basic 
equation (2.13) and the required boundary condition at z — 0, Eq.(3.4). 

(2) 

Asymptotic behavior of at z = oo is now easy to evaluate. Again, we find 
the behavior in agreement with general discussion of §3: 



qip In z — (\nz) 2 , 



p(2) _ _(2) 
r 2 — P2 



Qf = ,f-p«ln,, 



(4.14) 



where p^p an d can be read off from the general formula (4.5) in the previous 
subsection, 

5 



7T 

2 ' 



- 7 -ln2, 



(2) (2) 

and the constants p 2 an d % are found to be 



(2) 1 / m\ 2 1 /5\ 5 o 107, , 
F2 2 V 2 J 2\3J 24 210 u ; 



(2) 
9 2 



7T 

2 



?2 



107 
210 



(4.15) 



(4.16) 



Hence A| n to 0(e 2 ) is found to be 



a in 
A 2 



"2 ?e 



-ie In e 



,7T 



— 2- 



25 


5 






+ 18 + 


24 


107X1 




210^/ 


+ 



107 
210 



(4.17) 



5 , 107 

exp I ?e ?e 7r 

' 1 3 420 



(, n 2 r 25 5 2 io7 , , ON i 

( 1 - £ 2 +e {l8 + 24^ + 210^ + k2 )} + ' 
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where, in the last expression, we have absorbed the terms arising from qip and 
into the phase. 

(b) £ = 3 



f(1) 
'3 



,(1) 
'3 



(4.18) 



(2) 

Now we turn to the case of £ = 3. We focus on the real part of £3 . Then 
Eq.(2.16) yields 

1 A 1 ) , • f A 12 5 U( 

=-/3 + J3 / dz{2zn 3 - —n A - -n 2 )f$ 

/12 5 

+ 33 J dzn 3 j 3 -n 3 J dzj\ - ^' 3 (lnz) 2 . 

As before, after a lengthy calculation, we obtain 

c(2 ) 1 323 5065 /1031 445 \ 65 65 
^ =Tz M + Wz J2 -^ Jl -[^z + T^) JO + 6^ n °-6-z ni 

N / 3 13 9 15 13 . 
+ C(z) (^-n 3 - -n 2 - - ni - - - j3 

„, , / 3 . 13 . 9 . 15 . 13 
- S(z) ^-J 3 - y j2 - Y n - ^JO + T2 n, 

13 1 

+ B n (z)j 3 - Bj(z)n s - — j 3 In 2; - -j 3 (\nz) 2 + a^j 3 , 

(2) 

where a 3 is an arbitrary constant which will be put to zero hereafer. We have 
also examined the above result with Mathematica and confirmed it. 

(2) 

Asymptotic behavior of f 3 at z — 00 is 

fP^P^js + Q^ns; 

P^-pf + ihnz-ilnzf, ( 42 °) 



(4.19) 
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where and are found to be 



(2) 

Ps = 






(2) 
% = 


7T 

~2 





2 



2 1 /13V 5 2 13 



, H 7TH (7 + In 2), 

2 \ 6 / 24 42 u ; ' 



(4.21) 
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with and g^, again from Eq.(4.5), being given by 



rf^-f. ^ = f-7-ln2. (422) 



By absorbing the and 9^ terms into the phase, A| n to 0(e 2 ) is found to be 



A in J_ e -Hl,2e +1 ) ( . £ 13 _ 2 13 \ 

6 2 V 6 84 / 



x 



/ 7T 9 fl69 5 o 13. , ,1 



(4.23) 



(c) £ = 4 



The calculation of Q is parallel to the previous two cases. From Eq.(2.16), 
we have 

/f'=Re(«f) 

=-/i 1) + .74 / dz(2zn A - ln 5 - \n^)f^ 
Z J 6 6 

f 7 4m (424) 

- ra 4 / ^(22j4 - -J5 - gj'3) 7 4 

+ j 4 y dzn^-n^ y cbjf - ij 4 (lnz) 2 . 
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Then we obtain 



(2) _ 56 . / 5036 30334 V / 35252 _ 30334 14401 \ 
Ia ~165^ J5 + V 33V + 1155^2; J4+ V33^ _ 165V + 3465~I J 33 

/5036 45461 36287 \ . /140 5 \ 49 
~ VlT^ + 693i3 + 9240^y' Jl+ \^ ~ I8~z ) U ° ~ 6^ 

„, V 21 149 5 105 105 10 1571 N 

\10z 60 3z 2^ z d z 6930 , 

. . Z' 21 149 5 105 105 10 1571 \ 
- S{Z) Wz H - W + Yz n ~ 2V Jl - + T J0 + 6930™ 



1571 1 
+ B n (z) j 4 - Bj(z) n 4 - j a lnz ~2 j4 ( ln z ^ + , 

(2) 

where, for convenience, we shall choose a\ = also in this case. 
From the asymptotic behavior of at z = oo, we find 



(4.25) 



(2) 1 ( (i)\ 2 1 /149\ 2 5 2 1571 . , o . 
^ = -2^) + 2{-m) + 24^ + 6930 (7 + ln2) ' 

(4.26) 



2 V 6930 J ' 
where p 4 "^ and are once again read off from Eq.(4.5) as 



„?> = -!, ? « = ^- 7 -l„2. (4.27) 



Hence A 4 n to 0(e 2 ) is obtained as 



Af J-ie-^^) exp ( *e— - l e 2 —n) 
4 2 P V 60 13860 J 

(4.28) 

/ 7T 2 f 22201 5 2 1571, , \ 
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5. Conclusion 



We have formulated a new post-Newtonian type expansion to solve the ho- 
mogeneous Regge- Wheeler equation. We have then applied it to solve for the 
ingoing- wave Regge- Wheeler function Xf 1 , which is purely ingoing at the black- 
hole horizon. We have focused on X^ n since it plays an essential role in evaluating 
the gravitational waveform and luminosity radiated to infinity by a particle or- 
biting around a non-rotating black hole. The method has an advantage that the 
ingoing-wave boundary condition of Xj n at the horizon is naturally taken into ac- 
count. It also allows a systematic expansion of X\ n in powers of e = 2Muj, which 
is of 0(f 3 ) where v is the orbital velocity of a particle. 

By investigating general features of X^ n , we have argued that as long as we 
are concerned with the radiation going out to infinity, the effect of the radiation 
absorbed into the black-hole horizon is negligible until we go to 0(e 6 ), or 0(v 18 ) 
beyond Newtonian. Hence Xf 1 can be also used for a situation in which there is a 
non-rotating compact star instead of a black hole. Thus our method can provide 
an important guideline for the study of inspiraling compact binaries, though it is 
restricted to the case when one body has a small mass relative to the other. 

Then we have solved for Xf 1 and derived its analytic expression correct to 
0(e) for arbitrary I and to 0(e 2 ) for i = 2, 3 and 4. These results can be used to 
analytically evaluate the gravitational waveform and the luminosity radiated by a 
particle in circular orbit to 0(v 8 ) beyond Newtonian. This is in progress and the 

[15] 

result will be published in a forthcoming paper. 
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Appendix 

In this Appendix, we summarize mathematical formulas associated with the 
spherical Bessel functions which are used in the calculation of the ingoing-wave 
function XJ 1 . In doing so, we also sketch the calculation to 0(e 2 ). We denote any 
of the spherical Bessel functions by ( m and/or by 

The basic recurrence relation is 

2m + 1 

Cm-l + Cm+1 = Cm , (A.lJ 

Z 

from which it follows that 

i_ m = (-l) m n m _i <- n_ m = (-I)™" 1 j m -i • (A.2) 

Thus all the spherical Bessel functions appearing in the text can be expressed in 
terms of j m alone if necessary. There is also a useful series expansion formula: 

\2J ^— ^ r! (n — r)\ T{m + ra + r + 3/2) 

This enables us to get rid of all the inverse powers of z appearing, e.g., in the 
integrands of Eq.(4.12). Hence except for the terms involving Ci2z, Si2z and lnz, 
all the terms in the integrands reduce to either of the form zCmCn or of the form 
CmCn- For such forms, the following indefinite integral formulas are known to hold: 

/ dzCmCi = 7 w^— j 7T (CmC+1 + (m-l(n) CmC (m^n), 

J (m — n)(m + n + 1) V ,l+± / m — n 

J dzzCmCn = J dzzCm-lCn-1 ~ m + n (Cm-lCn-1 + CmCn) , 



dz Q( e 



2i 



j dzzQCi = j dzzCoC ~ y jcoCo + + mTl) ^ m + ' 



(A.4) 
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and 



(A.5) 



J dznl = - J dzjl - ^ = -S{z) - zr%, 

J dzn jo = -C(z) - zn jo - Inz , 

J dz zuq = — J dz zjq + In z = ~jC{z) + In z , 

J dzznojo = ~S{z), 
where the functions S(z) and C(z) are defined in Eq.(4.8), which we recapitulate: 
S(z) = Si2z, C(z) = Ci2z- 7 -m2z. (A.6) 

As for the terms involving C(z), S(z) and Inz, what we need to evaluate are 
the asymptotic forms of the integrals in the expressions of B n (z) and Bj(z) given 
by Eq.(4.9): 

2z 

/dx 
— sin x In x , 
x 





2z 



/dx 
— (1 - cosx) lnx, (A.7) 
x 





2z 

/dx 
— (Si x cos x — Ci x sin x) 

o 

First, let us evaluate I s . It can be done as 

2z 

dxx a ~ x sinx 



9 / , sinfa\ (A.8) 

r(i + «)- 2 



7T 



= "2 T 
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In a similar manner, I c is evaluated as 



d 



2z 



Ic.(z) = TT- ( / dxx a 1 (1 — cos x) 
da 



a=0 



d_ I {2z) a - 1 
da \ a 



2z 



I 



dxx a 1 cos x — 



a 



a=0 



(A.9) 



— ► -(\n2zY - — 

z^oo 2 v ; da 



8 /r(l + a)cosfa-l 



a 



a=0 



Finally, we consider I cs . We utilize the following integral representation of the 
integrand: 



si x cos x-cixsini = 



dt 



sint 
7+x' 



(A.10) 



where six = Si x — n/2 and cix = Cix. Then we manipulate the integral as 
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2z ( oo 

, , . , , sin t ( 1 1 \ 7r cos x 
hs= I dx{ / dt— ( — ) + -- 



KO 



t \ t + x x I 2 x 



oo 2z 2z 

dt . f dx 7i f , cos x — 1 



/ — sm t / h - / 

7 * J t+x 2J 





dx- 



X 



oo oo 

= J y sintln(2^ + t) - J jsintlnt+ ^C(z) (A.ll) 
o o 

oo 

= f (ln2z + 7 + C(z)) + | ysintln (l + ^) 



oo 

— sin2zt ln(l + t) — > 0. 

t 2^00 



Hence J cs vanishes as z — > oo. 
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